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ABSTRACT 

Many quasi-simultaneous optical observations of 17 blazars were obtained 
from the previous published papers over last 19 years in order to investi- 
gate the spectral slope variability and understand the radiation mechanism of 

o 

\^ • blazars. The long period dereddened optical spectral slopes were calculated 



in this paper. Our analysis upon the average spectral slope distribution sug- 
gests that the spectra of Flat Spectrum Radio Quasars (FSRQs) and High 

. energy peaked BL Lac objects (HBLs) are probably deformed by other emis- 

5h ' 

' sion components. The average spectral slopes of Low energy peaked BL Lac 

objects(LBLs), which scatter around 1.5, show a good accordance with Syn- 
chrotron Self-Compton (SSC) loss-dominated model. We present and discuss 
the variability between spectral slope and optical luminosity. The spectra of 
all HBLs and LBLs get flatter when they turn brighter, while for FSRQs this 
trend does not exist or may be in a reverse situation. This phenomenon may 
imply that there is a thermal contribution to the optical spectrum for FSRQs. 
For the FSRQ 1156+295, there is a hint that the slope gets flatter at both the 
brightest and faintest states. Our result shows that three subclasses locate in 
different regions in the pattern of slope variability indicator versus the average 
spectral slope. Furthermore we also delineate that relativistic jet mechanism 
is supported by the significant correlation between the optical Doppler factor 
and the average spectral slope. 



Key words: galaxies: active - BL Lacertae objects: general - quasars: general 
— galaxies: fundamental parameters. 
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1 INTRODUCTION 

Blazars are a subset of Active Galactic Nuclei (AGNs), which are compact, fiat spectrum ra- 
dio sources with highly variable and polarized nonther mal continuum emission ranging from 
radio up to X-ray and often to 7-ray frequency (e.g. 



Sillanpaa et al. 



1991 



Maraschi et al. 



199 



Fugmann 



Urrv fc Padovani 



1988; 



] 199.4 



Bregman et al 



1990; 



Their observed proper- 



ties a re mainly considered to originate from a relativistic beaming jet (jRees. Brgelman fc Blanford 
19811 ). which is possibly powered and accelerated by a rotating and accreting supermassive 
black hole. There are two classes of blazars, BL Lac objects (BL Lacs) and FSRQs, the 
former one has a featureless optical continuum, while the latter has a great many strong and 
broad emission lines. According to the peak frequency of synchrotron emission relative to the 
peak flux, BL Lacs are classified as LBLs, whose peak is located in infr ared or optical band. 



and 



5BLs, whose peak is l ocated in the ultraviolet or X-ray frequency (jPadovani fc Giommi 



Giommi 



1999ll. For HBL s. the b roadband spectral index anx is le ss than or 



equal to 0.75 (jUrrv fc Padovanilll995| ) or 0.80 (ISambruna. Maraschi fc Urrvlll996[ ). while for 
LBLs, it is greater than 0.75 or 0.80. 

The studies of the spectral energy distribution (SED) of blazars are very important for the 
understanding of the physical radiation mechanism and to constrain the parameters involved 
in physical model, and the studies are crucial in analyzing individual emission components. 
Synchrotron inverse Compton emission model predicts that the spectrum gets harder as the 
source turns brighter. Many investigators found that the variation amplitude of blazars at 
high frequencies was larger than that at low frequencies, namely, the spectrum became flatter 
when the flux increased, while the spectrum became steeper wh en the flux decreased fe.g . 



R.acine 



mi 



Gear. Robson fc Brown 



mm 



Massano 



sugges ted that it might not be always correct (e.g. 



e1^ 



Ramirez et al. 



19971). But iGhosh et al 



20041) . 



D'Amicis et al 



()2002| ) reported the optical spectral slope variability of 8 BL Lacs. They found that the 



Jiects. Four objects 



Trevese fc Vagnetti 



spectral slope became flatter when the source turned brighter for all 8 o 
showed marked correlation between spectral slopes and R magnitudes. 

analyzed the spectral slope variability of 42 PG quasars. They concluded that the 
spectral variability must be intrinsic of the nuclear component. Their numerical simulation 
showed that hot spots on the disk were able to represent the observed spectral variability, 
while single host galaxy contribution or single changes of accretion rate was insufficient to 
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D'Amicis et al. 



explai n the slope variability. IVagnetti et al.l ()2003l ) analyzed the data from 
( 20021 ) under a simple synchrotron mode and the thermal bump (TB) mode. They showed 
that the spectral variability, even restricted to the optical band, could be used to set limits 
on the relative contribution of the synchrotron component and the thermal component. They 
deduced that the spectral variability of BL Lacs differed from the variability of quasars. The 



8 BL Lacs they studied were all LBLs. 



Ramirez et al. 



( 2004^ presented that the spectrum 



of a FSRQ PKS 0736+017 became softer when it turned brighter. The correlation between 
spectral slope and flux was strong. Its variability relation between slope and flux was reverse 
to the result that was mentioned above for LBLs. 

In this paper we show and discuss the optical spectral slope variability of 17 blazars, 
which include 9 LBLs, 4 HBLs and 4 FSRQs. In section 2 we describe the data and data 
reduction process. The average optical spectral slope distribution is analyzed in section 3, 
the spectral slope variability is given in section 4, and conclusions are presented in the last 
section. 



2 DATA AND DATA REDUCTION 



We collected optical observations from the published papers in period of 1985 ~ 2003, which 
were made in the standard Johnson- Cousins B, V, R, I bands. The observations should be 
simultaneous in all above-mentioned optical bands i n order to get th e optical spectral slope. 



XieetaL 



(20051) reported that the min- 



but it is difficult due to telescopes and instruments, 
imum variability timescales of 21 blazars were longer than 1 hour, and most of them were 
longer than 2 or 3 hours. We regard this quasi-s imultaneous mul t i-ban d obs ervations taken 



and 



Fiorucci et al 



within 1 hour as "simultaneous" observations (jP'Amicis et al 

12004! ) accepted quasi- simultaneous observations taken within 2 hours). So we merely se- 
lected the quasi-simultaneous observations, which were made in three(any three bands of 
B, V, R and /) or four bands(-B, V, R, I) within 1 hour, as the raw data for accuracy. We 
selected a group of quasi-simultaneous observations in three or four bands as one data group, 
and most of the data groups were made within 30 minutes. Then the following reduction 
processes were applied to all the selected raw data. All the observation data made in V band 
were corrected for the foreground Galactic interstellar reddening and absorption with the 
employment of the e xtinction value Ay, which is d e duced from the maps of dust infrared 
emission reported bv ISchlegel. Finkbeiner fc David (|l998h . The extinction values in other 
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bands were obtained by using the curve of ICardelli. Clavton fc Mathid (|l989| ). The red- 
dening and absorption correction can be done by subtracting the corresponding extinction 
value from the raw data. The conversion from apparent magnitudes to optical flux densitie s 
was made by using the zero magnitude equivalent flux density given by 



Mead et al. 



(19901). 



The rest frequenci es were obtained through the effective wavelengths of t he filters given by 



Mead et al. 



( 19901 1 and the redshift given by IVeron-Cettv fc VeronI (|2003[ ) . 



The radiation of blazars can be described by a single power law = Av~" in optical 
band, where a is the spectral slope. It was obtained by using the linear least square fitting 
between logF^^ and logv. We can obtain one a from each data group, but we accepted only 
the spectral slope whose square of Pearson's linear correlation coefficient is greater than 0.9 
and whose standard deviation of fitting slope is less than 0.4 (the typical standard deviation 
of the spectral slope is less than 0.1). The aforesaid selection criterions assure the accuracy of 



the sl opes by cutting off the probably wrong values due to quasi-simultaneous (|Fiorucci et al 



2004f ). After all the above reduction, we accepted 1418 useful data groups from all the raw 



data groups for 17 blazar samples whose data were sufficient to investigate the spectral slope 
variability. Our main results are li sted in Tabled The sequence columns in Table d are the 
source name; the redshift z from IVeron-Cettv fc VeronI ()2003r i: the average R magnitude 
R and the variable amplitude (max-min) of R magnitude Mr\ the average optical spectral 
slope Uave and the standard deviation of aave'-, the variable amplitude(max-min) of the optical 
spectral slope Ma-, the slope of the linear regression between a and the R magnitudes, 6, (it 
is taken as the slope variability indicator), followed by its standard deviation; the correlation 
coefficient between a and R magnitude, r^; the optical Doppler factor 5; the reference of 5; 
the extinction value in V band Ay; subclass of the source; the number of useful data groups; 
the references of optical photometric observations. 



3 AVERAGE OPTICAL SPECTRAL SLOPE DISTRIBUTION 



Blazars are classified into three subclasses (Donato et al 



2001 



Padovani fc Giommi 



19951). 



This classification is based upon the physical natures and characters. So physical connections 
among these three subclasses should exist in manv fields. Manv investigators studied such 



connections (iGiommi et 



1996 



Hi .11 9901 : 



D'Eha fc Cavahere 



Lamer. Brunner fc Staubertll99fi 



Stecker. de .Tager Salamon 



20011 ). which would play an important role in constraini ng the radi- 



ation mechanism of blazars and understanding the fundamental nature of blazars. 



Fossati et al. 
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Table 1. Principal data and results of 17 blazar samples 
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( 19981 ) showed that blazars form a sequence, the SED changed continuously with the change 
of their bolometric luminosity. The synchrotron component and inverse Compton compo- 
nent of HBLs approximately had the same scale of energy and showed two emission peaks, 
which shifted to lower frequencies as the bolometric luminosity increases. Subclasses from 
FSRQs through LBLs to HBLs exhibited a sequenc e of increasing spectral hardness with the 
decreasing of luminosity ( Bottcher k. Dermeiil2002t ). It suggests that this may be a sequence 
of decreasing accretion rate along the same sequence, so they proposed an evolutionary sce- 
nario linking FSRQs, LBLs and HBLs through the gradual depletion of the circumnuclear 
environment of a supermassive black hole. This trend provides the opportunity to unite 
three subclasses of blazars under one single scheme. 

Analyzing the dereddened optical spectrum can distinguish the optical radiation com- 
ponent of three subclasses of blazars properly. The average optical s pectral slopes are listed 



in cloumn 4 of Tabled Our results are consistent with the results of 



D'Amicis et al. 



(2002) 



for 6 commo n LBLs within 
the results of iFiorucci et al 



(jfst andard deviation), and they are in good agreement with 
()2004i l for 14 common blazars within 1 a, but the adopted pho- 



tometric observations are different. The distribution of aave for our 17 samples is given in 
Fig. [T] by subclass. The average of aave is 1.07, 1.49 and 1.08 for FSRQs, LBLs and HBLs, 
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Figure 1. Histogram of the average optical spectral slopes for our 17 samples by subclass 



respectively. We give the distribution of an,,^ in Fig 



i nclud e our 17 samples, 8 samples frorn 



( 20041 ). 37 objects from 



Fiorucci et a1. 



3 for 4 1 blazars by s ubclass, which 



D'Amicis et al. 



( 20021 ). 1 FSRQ from 



Ramirez et al. 



I2OO4 ). We averaged the spectral slopes for all the 



same available sources. The average of aave is 1.10, 1.60 and 1.18 for FSRQs, LBLs and 
HBLs, respectively. The spectra of LBLs are obviously steeper than those of HBLs and 
FSRQs. According to the peak frequency of synchrotron emission, we expect aave > 1 for 
LBLs and FSRQs due to their SED in the descending part, b ut we expect an„p < 1 for 



HBLs because of their SED in the ascending part. The model of iChiang &: Bottcheii (j2002| ) 
predicted that for a broad range of particle injection distribution. Synchrotron Self Comp- 
ton (SSC)-loss-dominated synchrotron emission exhibited a spectrum whose spectral slope 
was 1.5 in optical band. The histogram in Fig. |21 shows that all aave of LBLs are greater 
than 1 and the distribution peak is around 1.5 (average is 1.60). Our result is consistent 
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Figure 2. Histogram of tlie average optical spectral slopes for 41 samples by subclass 



with 



Fiorucci et al 



( 20041 ) ■ it indicates that SSC- loss-dominated mechanism is probably the 
main radiation mechanism for LBLs. aave of HBLs distribute in the range from 0.5 to 2, its 
distribution peak is between 0.75 and 1.25, while the predicted spectral slope under pure 
synchrotron emission should be less than or equal to 1. This means that the optical radi- 
ation of HBLs may contain other components, such as thermal contribution of accretion 
disk or host galaxy, nonthermal emission originating from different regions or different sizes 
of the relativistic jet. For FSRQs, the distribution range of aave is 0.5 ~ 2, whose peak is 
between 0.75 ~ 1, so their optical emission is likely to be deformed by thermal emission 
or other components. Synchrotron component does not domina te and other mechan ism can 



Turler et al 



( 1999[ l showed 



not be negligible. "Thermal Bump" (TB) may be an example, 
the "UV bump" for 3C 273. More simultaneous multi-band observations of larger samples 
are important to thoroughly study the emission components of blazars. 
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Figure 3. The optical spectral slope versus R magnitude for 9 LBLs 



4 SPECTRAL SLOPE VARIABILITY 

The spectral slopes and their uncertainties were obtained by linear least square fitting 
with the "simultaneous" raw data groups. Figs. show the variability relations be- 



tween the spectral slopes and the R magnitudes (jMassaro fc Trevesd Il996a[ ) for 9 LBLs, 
4 HBLs and 4 FSRQs, respectively. Linear regression analysis was applied to every sam- 
ple in order to investigate the possible relation between the optical spectral slope and the 
source luminosity. The linear regression slopes b, their standard deviations and the cor- 
relation coefficients are listed in Table ^ A tendency to become ffat with increasing 
of R magnitude can be observed for all LBLs and HBLs. are greater than 0.5 for 3 
out of 4 HBLs and for 5 out of 9 LBLs. The average of correlation coefficient is 0.60 and 
0.54 for HBLs and LBLs, respectively. The correlation relations are significant at a 0.95 
confidence level for all HBLs and LBLs except for 1402+042, the average correlation coef- 
ficient of HBLs is higher than that of LBLs. The variability is consistent with the relation 
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Figure 4. The optical spectral slope versus R magnitude for 4 HBLs 



variability relation in other wavebands (e.g. Gear. Robson & Brown 1986 


Massaro et al 




199S; 


Kedziora-Chudczer et aLlLoOl: 


Panadakis et ahlLoO 


a) . Many models ( C 


ian et al .111 991 


T.i <^ KusimosJbnnn: 


Snada et alhoni : Sikora et al. 2001 


Wane- & Kusunosf 


I2OO2: Chiang & Rottcher 


2002; 


Ivaffnetti et al. 


2003) stated that spectra of blazars became flat when they turned 



bright. For FSRQs, the spectra generally become 



1981 



Brown et al. 



1986 



steep 



Clements. Jenks fc Torres 



2003 



when they turn brig ht (e.g. 



Ramirez et al. 



Miller 



20041 ) except for 3C 



279, but all four correlations are not significant at a 0.95 confidence level. It is interesting 
that the tendency is possibly opposite in the faint state against the bri ght state for 1156+295 . 
The general result can be interpreted by the thermal bump model bv lVagnetti et al.l ()2003r i. 
they predicted that a strong thermal bump could produce a spectrum that becomes steep 
when getting bright, where the thermal bump emerges in the faint state. The variability of 
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Figure 5. The optical spectral slope versus R magnitude for 4 FSRQs 



1156+295 is a good support to 



Brown et al, 



they mentioned that in high bright- 



ness state the spectra variabihty of FSRQs showed similar properties with that of BL Lacs. 
But in the faint state, since the thermal contribution was larger in the blue part, the nonther- 
mal component was steeper than the composite spectrum. The nonthermal component had 
an even more dominant contrib ution to the total flu x when the object got bright, and the 



composite spectrum got steeper (jRamirez et al 



2004 ). Our results suggest that the emission 



of FSRQs probably contains thermal component and nonthermal synchrotron component, 
and the former one flattens the spectrum. 

The relation between the regression slope b (it is the optical spectral slope variability 
indictor) and aave is shown in Fig. IHl We can see that different subclasses have the tendency 
to locate in different regions in this pattern, LBLs are between FSRQs and HBLs . This 
may connect with the evolutionary scenario obtained bv iBottcher fc Dermerl ()2002i l. The 
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accretion rate decreases along the sequence of FSRQs, LBLs and HBLs, due to the gradual 
depletion of the circumnuclear environment of a supermassive black hole. More samples are 
needed to confirm this result. The mean value of b is —0.24, 0.32 and 1.12 for FSRQs, 
LBLs and HBLs, respectively. The slope variability indicator of HBLs is the greatest, but 
the average variable amplitude of R magnitudes is the smallest and the average variable 
amplitude of the spectral slope is smaller than that of LBLs. 



The observed flux density is enhanced by relativistic Doppler beaming effect. 



Ghisellini &: Maraschi 



119891 ) presented that the jet plasma might accelerate outward, so that the bulk velocity 
and the bulk Lo rentz factor incr eased with incr e asing distance along the jet. According 



Fan et al 



( 19931 ) and IXie et al.l (|200lh obtained an empirical formula of 



to this proposal, 

frequency-dependent Doppler factor, 6^, = 5^+^/^^°9{'^o/'^) ^ where u is the frequency of con- 
sidered waveband, pq is the frequency of optical band, 5 is the optical Doppler factor, 5y 
is the Doppler factor at frequency v. Therefore the optical Doppler factor can be obtained 
if the radio or X-ray Doppler factor is known. Fig. [3 shows the relation between aave and 
the optical Doppler factor 5. The linear regression analysis between 5 and aave shows that 
the linear correlation is significant at a 0.99 confidence level. The correlation coefficient is 
0.79, and the regression line is plotted with solid fine. In order to exclude the dependence 
on the outlying point 0235-1-164, we applied linear regression analysis to the data excluded 
0235-1-164. The regression coefficient is consistent with the above-mentioned results within 
one standard deviation, the correlation coefficient is 0.70, the linear correlation is also sig- 
nificant at a 0.99 confidence level. The regression line is shown with dotted-line in Fig. [3 
Two regression lines are almost overlapped, thus the result is not dependent on the outlying 
point significantly. The regression equation for all the data is: 

aave = 0.573(±0.084) + 0.297(±0.045)5 (1) 

We can obtain the optical Doppler factor from average optical spectral slope if this relation 
is confirmed. This result suggests that Doppler beaming effect of relativistic jet is the main 
mechanism for the understanding of the properties of blazars. 



5 CONCLUSIONS 

We collected a large number of "simultaneous" {At < 1 hour and most of their At < 30 
minutes) optical photometric data in a long period from 1985 to 2003. Here we gathered a 
larger blazar sample which includes some objects from each of three subclasses. A compar- 
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Figure 6. The average optical spectral slope aave versus the slope variability indicator (linear regression slope) b 

ison between these subclasses enables us to gain some useful information on the emission 
mechanisms and emission components. The research of this sample is a great benefit to 
investigate the natures of blazars. In this paper we calculated the long period dereddened 
optical spectral slopes by linear least square fitting, which is based on the power law relation. 
The average optical spectral slope distributions of three subclasses of blazars were analyzed. 
We further on analyzed the relation between the optical spectral slope and R magnitude, 
the relation between average optical spectral slope and slope variability indicator b, and 
the relation between average optical spectral slope and optical Doppler factor. A discussion 
mainly about emission mechanism, emission components and connection of three subclasses 
of blazars presents the following conclusions: 

(i) All the average spectral slopes of our 9 LBLs are greater than unity and scatter aro und 



1.5, this is a good support to the SSC-loss-dominated model ([Chiang fc Bottcheill2002l ). 
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Figure 7. The relation between average optical spectral slope aave and optical Doppler factor S. Solid line is the linear 
regression to all the data, and dotted-line is the linear regression to the data except for the outlying point 0235+164. 

(ii) The optical spectra of all our HBLs and LBLs become flatter when the sources turn 
brighter, the average linear regression correlation coefficient of HBLs is higher than that of 
LBLs. The relation between spectral slope and R magnitude of FSRQs suggests that the 
emission of FSRQs probably has "thermal bump" contribution and the thermal component 
must be considered. 

(iii) The slope distribution and variability for HBLs indicate that their spectra appear to 
be deformed by other components, which presumably come from the thermal accretion disk 
or originate from different regions of the jet. 

(iv) Different subclasses of blazars have the tendency to locate in different regions in the 
pattern of the slope variability indicator versus average spectral slope. 

(v) There is a significant correlation (correlation coefficient is 0.79) between the average 
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optical spectral slope and optical Doppler factor, which suggests that Doppler beaming effect 
of relativistic jet is the main mechanism for blazars. 
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